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INTERFERENCE REDUCTION BY STEP FUNCTION 

REMOVAL 

TECHNICAL FIELD 

This invention relates to reception of a signal. 

BACKGROUND 

A received signal may include a desired signal from a 
desired source along with one or more undesired signals, such 
as, for example, a noise signal (e.g., additive noise such as 
white Gaussian noise) from a noise source, and/or an 
interfering signal from an interfering source (e.g., main-lobe 
or side-lobe energy of the interfering signal) . The received 
signal also may include an offset component such as, for 
example, a DC (direct current) offset component, that may be 
undesirable. The offset component is an additional additive 
term and may be a constant offset such as, for example, a DC 
offset, or may be a non-constant offset such as, for example, 
a step function. 

To extract the desired signal from the received signal, 

characteristic parameters (e.g., data bits, frequency offset, 

DC offset) that model the received signal may be estimated. 

It may be desirable to perform preprocessing of the received 

signal prior to estimating the characteristic parameters, such 
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as, for example, estimating the offset and removing its 
effect. For example, the offset may be estimated as a mean of 
the received signal and the mean may then be subtracted from 
the received signal . 

The offset of the received signal may vary, for example, 
because of variation in the interfering signal. Such 
variations may cause the mean of the received signal to 
provide a poor estimate of the signal offset. Subtracting a 
poor estimate of the offset would then bias the estimates of 
the characteristic parameters and lead to inaccurate results. 

DESCRIPTION OF DRAWINGS 

Fig. 1 is a schematic diagram of a communication system 
configured to estimate and correct a signal having an offset 
that may be modeled as a step function. 

Fig. 2 is a diagram illustrating two time-division, 
multiple-access (TDMA) users that are not aligned in time and 
that may cause interference which may appear as an additional 
signal offset in the form of a step function. 

Fig. 3 is a schematic diagram of a receiver that may be 
used with the communication system of Fig. 1. 

Fig. 4 is a diagram illustrating a step function that may 
be used to model the signal offset of the system of Fig. 1. 



Fig. 5 is a schematic flow diagram illustrating a 
systematic process for offset correcting a signal having an 
offset that may be modeled as a step function by removing the 
step function of Fig. 4 from the received signal. 

Like reference symbols in the various drawings indicate 
like elements. 

DETAILED DESCRIPTION 

For illustrative purposes, a process is described for 
interference reduction by offset correcting a signal, where 
the offset may be modeled as including a step function and the 
signal is corrected by removing the undesired step function. 
For clarity of exposition, the description generally proceeds 
from an account of general elements and their high level 
relationship to a detailed account of illustrative roles, 
configurations, and components of the elements. 

Referring to Fig. 1, a generalized system 100 (e.g., a 
global system for mobile communications (GSM) , a time- 
division, multiple-access (TDMA) system, or a frequency- 
division, multiple access (FDMA) system) may be used to 
receive a transmitted signal and to correct the offset of the 
received signal, where the signal offset may be modeled as 
including a step function. Exemplary components of the system 
100 are described in greater detail below. 



The system 100 of Fig. 1 generally includes a transmitter 
110, a receiver 130 (e.g., a superheterodyne receiver, a dual- 
conversion superheterodyne receiver, or a direct conversion 
receiver) , and a channel 150 that models how the environment 
has changed the transmitted signal as perceived at the input 
port of the receiver. 

In general, the transmitter 110 and the receiver 13 0 may 
include any devices, systems, pieces of code, and/or 
combinations of these that may be used to transmit or receive, 
respectively, a waveform z(t) that generally may be 
represented as 

z(0 = Re{s(t)}cos(co 0 t) - \m{s(t)}sm(co Q t) . (l . o) 

In equation (1.0), s(t) may denote a complex signal, co 0 = 2nf 0 
may be an associated carrier frequency, and Re{s(t)} and 
Im{s(t)} denote respectively, real and imaginary parts of 
s(t) . 

A transmitter (e.g., transmitter 110) and/or a receiver 
(e.g., receiver 130) generally may include, for example, a 
mixer (e.g., mixer 135), a summer, a phase locked loop, a 
frequency synthesizer, a filter (e.g., filter 137), an 
oscillator (e.g., local oscillator 133), a frequency divider, 
a phase modulator, a down converter, an amplifier, a phase 
shifter, an analog-to-digital (A/D) converter or a digital-to- 
analog (D/A) converter (e.g., A/D converter 137), a 



microprocessor (MPU) , a digital signal processor (DSP) , a 
computer, or a signal processing circuit, whether linear or 
nonlinear, analog or digital, and/or any combination of these 
elements . 

More specifically, receiver 130 may include a down- 
converter for down -converting an input signal from radio 
frequency (RF) to baseband. The down- converter includes the 
local oscillator 133, the mixer 135, and the filter 137 (e.g., 
an infinite impulse response filter, a finite impulse response 
filter) . The receiver 130 also may include an A/D converter 
to generate a discrete signal from a continuous input and, for 
performing step parameter estimation and offset correction, 
any device, system, or piece of code suitable for that task, 
such as, for example, step parameter estimator and signal 
offset corrector 13 9. The step parameter estimator and signal 
offset corrector 13 9 may include, for example, a 
microprocessor control unit (MCU) , a digital signal processing 
(DSP) component, a computer, a piece of code, a signal 
processing circuit, whether linear or nonlinear, analog or 
digital, and/or any combination of these for use in performing 
the step parameter estimation and/or the offset correction, 
including step function removal. 

The transmitter 110 transmits a signal z(t) over the 
channel 150. The channel 150 may include any medium over 



which a signal may be communicated, such as, for example, an 
RF (radio frequency) portion of the electromagnetic spectrum, 
and or any other portion of the electromagnetic spectrum. 
Associated with the channel are an interference source 151 
that generates an interference signal J(t) and a noise source 
153 that generates a noise signal w(t) . The noise source 153 
and the interference source 151 add noise w(t) and 
interference J(t), respectively, to z(t) to form a signal r(t) 
received by the receiver. 

The noise w(t) may include, for example, additive white 
Gaussian noise that may have a zero or non-zero mean, while 
the interfering signal J(t) may have very different 
characteristics before and after an event that occurs within 
the burst. For example, if the interfering source is due to a 
different TDMA user who is transmitting at the same frequency 
as, but not time aligned with, the desired user, then the 
interference may appear as being turned on and off during the 
burst for multiple bursts. 

Fig. 2 illustrates one example of burst interference that 
may generate a signal offset that may be represented as 
including a step function at the input of the data bit 
estimator 140. The data bit estimator 140 may include, for 
example, a matched filter, and/or a decoder such as, for 
example, a convolution decoder, and may perform functions 



including de-interleaving or decoding, and further may provide 
an estimate of data bits sent over the channel 150. 

Note that the signals z(t) and T(t), as shown in Fig. 2, 
are illustrative only and may not represent certain 
characteristics of actual physical signals. As shown, signal 
z(t) is transmitted in time slot n of a first TDMA waveform 
210, while the interference signal J(t) is transmitted in time 
slot m of a second TDMA waveform 220. Each TDMA waveform is 
associated with different TDMA channels (e.g., different GSM 
base stations with or without different hopping patterns) . 

The time slots for these first and second TDMA waveforms 
are not time aligned with each other (e.g., each time slot of 
the second TDMA waveform lags (or leads) the corresponding 
time slot of the first TDMA channel by the same time increment 
of t 2 -ti) . The interference signal J(t) also may have a power 
that is much greater than that of z(t) and a center frequency 
different than the center frequency co 0 of z(t), such as, for 
example, a center frequency that approximates a harmonic of 
co 0 . 

Because of the phase difference between the two TDMA 
waveforms, the transmission of signal J(t) at time t 2 may 
appear as interference that is turned on and off and is 
included in signal z(t). Moreover, a TDMA channel structure, 
such as, for example, a TDMA time-slot assignment methodology, 




may ensure that z(t) and I(t) transmit in lockstep, causing 
z(t) to experience burst interference from J(t) beginning at 
the same relative point in each time-slot in which z(t) is 
transmitted (e.g., t 2 -ti from the beginning of each time slot) . 
5 Referring again to Fig. 1, the receiver 13 0 receives from 

the channel a signal r(t) that includes z(t) + I(t) + w(t) . A 
mixer 135 produces y(t) by mixing r(t) with the combination of 
a sinusoidal signal Q{t) generated by local oscillator 133 

l* with an attenuated version of r(t) that may leak into the 

f i 

O 10 local oscillator 133. The leakage of r(t) into the local 

■P 

y oscillator 133 is represented by multiplying the received 

^ signal r(t) by an attenuation factor y to produce yr(t), and 

jj\ then summing yr(t) with the output Q(t) of an ideal local 

rU 

flj oscillator 134. Leakage of r(t) into the local oscillator 133 

p 15 causes y(t) to include the signal mix of r(t) [Q{t) + yr(t)] . 

Thereafter, y(t) passes through a low pass or band pass 
filter and/or an A/D (analog-to-digital) converter 137 (e.g., 
an integrator that performs the functions of A/D conversion 
and low pass filtering) to produce a discrete signal y n that 
20 may include an undesirable offset component. Thereafter, y n 
is processed further by step parameter estimator and offset 
corrector 139, which models the offset as a step function and 
estimates parameters descriptive of the step function. The 
offset of y n is corrected by offset corrector 139 based on the 
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estimated step function parameters. 

Fig. 3 illustrates a receiver 130 that may be used to 
implement the system of Fig. 1, and in which a signal r(t) 
leaks into an ideal local oscillator 334. The signal r(t) may 
include a transmitted signal z(t) , a sum of interfering 
signals T(t), and additive white Gaussian noise signal w{t) . 
Signal z(t) may be represented as a real signal resulting from 
upconversion of a complex signal s(t): 

zft;=Re{s (t) }cos cD 0 t-Im{s (t) }sin co 0 t. (1.1) 

Due to the leakage of r(t) into the local oscillator 333, the 
mixer 335 may not simply multiply r(t) by a sinusoid (e .q . f A 0 e~ Ja)ot , 
where A 0 is a known value) . Instead r(t) is multiplied by the 
sinusoid and an attenuated version of the input, ^r(t) . The 
resulting signal may be expressed as: 

y(t) = A,[z(t) + 7(0 + w(0]e" M ' + yl\t) + yz\t) + yw\t) + 
2yz{t)I(t) + 2yz(t)w(t) + 2M0'(0- 

Substituting equation 1.1 for z(t) produces 
y{t) = 4[Re{5(0}cos<2> 0 ( - lm{s(0}sin# 0 /]e~ >0 ' 

+ Aj{t)e Jm * + 4M0e~ >0 ' + yi\t) 

+ /[Refc(0}cosay - ]m{s(t)}sma) 0 tf + yw 2 (t) (1.3) 

+ 2/[Re{ t s(0}coscy 0 ^ - Im{s(/)}sina> 0 /]/(0 

+ 2/[Re{s(0}costfV " bn{s(t)}smo) Q t]w(t) + 2yw(f)I{t). 

If the attenuation term y is sufficiently small compared to the 

signal amplitude, then y(t) may be approximated as 



# 

y(t) « ^[Re^COjcos^ - Im^COjsin^K 7 " 0 ' 

+ + 4M'K >0 ' + yi\t\ ( 1 ■ 4 ) 

in which the term y j2 ( t ) is retained because it is assumed that 
J(t) is of substantially greater power than z(t) . 

The signal y(t) then passes through a low pass filter 339, for 
example, to produce yi ow (t) , where yi ov (t) may be approximated as: 

y ibw (0 * y*(0 + ^ (0 + w w (0 , (i.5) 

in which the term A Q I(t)e' JO)ot is assumed to be substantially removed 
by the low pass filter and, therefore, has been dropped. In 
equation (1.5), the term w bb (t) represents a baseband portion of 
A Q w(t)e~ J6>ot that remains after passage through the low pass filter 

339. Assuming J(t) to be generally sinusoidal, yl 2 {t) may include an 
offset component (e.g., a DC offset) and a bandpass component at 
twice the center frequency of J(t) . The low pass filter may 
substantially remove the bandpass component of yl 2 (t) while leaving 
essentially unaffected the offset component, represented in 
equation (1.5) as yl\ b . When for example, I(t) is switched on or 
off, the offset component of y^lb ma Y be modeled as a step function. 

An A/D converter 341 may be used to generate a discrete 
signal y n based on the signal y low (t) . Assuming that yl 2 bb may 

be represented as a step function, the discrete signal y n may 
be represented as: 
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y n * -f s. (0) + c\ + (c2 - clK_ a + h>„ (1.6) 

where s n (0; is a discrete model of the baseband signal, 0 is a 
vector of unknown signal parameters (e.g., data bits, frequency- 
offset) , and w n is a discrete representation of zero-mean 
additive white Gaussian noise remaining after passing w(t) 
through the low pass filter 339 and the A/D converter 341. 
Also, referring now to Fig. 4, u n represents a unit step 
function that transitions from zero to one at n equals zero, 
such that cl + (c2 - c\)u n _ a represents a step function with amplitude 
of cl before the step transition and amplitude of c2 after the 
step transition, where the step transition occurs at time n 
equals a. 

Referring again to Fig. 3, the signal y n is provided to the 
parameter estimator and offset corrector 139. The parameters 
cl, c2 and a of the step function are estimated, and the 
estimated parameters then are used to correct the offset of 
signal y n , to produce an output signal that may be represented 
as : 

-f S »W + W >> ■ (1.7) 

The parameter estimator 13 9 may estimate the step function 
parameters based on, for example, gradient descent algorithms 
(e.g., the least-mean-square algorithm, Newton's method, the 
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steepest descent method, and/or any combination of these methods) 
and/or the maximum- likelihood (ML) method. 

The ML method provides a general method of maximizing the 
likelihood of the joint probability density function of the 
values of the received signal vector (y 0 , . . . f y N -i) given an 
intended signal vector (x 0 , . . . ,*n-i) . For the case when the 
observations are independent, a combined probability, or 
likelihood function, may be expressed as the product of the 
probability densities of the independent received signal vector 
samples, i.e., p = p(y 0 )^^p(y N _ l ) , where it may be assumed that 
each probability density can be modeled as a Gaussian density. 
The likelihood function p is then maximized to find the optimal 
parameters using any suitable optimization technique (e.g., a 
non-linear optimization technique) , such as, for example, the 
Nelder-Mead method (a method based upon the simplex algorithm) , 
the steepest descent method, the LMS (least-mean-square) 
method, the Levenberg-Marquardt method (a least squares 
approach) , the Davidson-Fletcher-Powell method (a quasi-Newton 
based method) , or the Broyden-Fletcher-Goldf arb-Shannon method 
(a quasi-Newton based method) , and/or any combination of one or 
more of these or other optimization methods. 

More specifically, a ML estimate of the step function 
parameters cl, c2 , and a can be obtained from the samples y n as 

described above. For example, we may take the baseband signal 
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model s n {9) and noise model w(n) to have a zero mean, since 
their means can be incorporated into the step function 

parameters. Using ^-5 n (^) + cl + (c2-cl)w n _ a as an expression of the 

mean of the individual values of the received signal vector 
produces the following ML likelihood function of the complex 
observation: 

£4 1 -\y„-X n (0^c\Hc2-c\)u n J L 2 

p=T\-r— e < (1 - 8 > 

n=0 y/7T<J 

which may be simplified to 

1 



P = 



f , \" _J_£L_A 5ii(<?)+c1+(c2 _ c1)i||i . 



(1.9) 



\\jncF J 

To maximize the value of p, it is sufficient to minimize 
the value of 

N-\ 



/-z 



n=0 



yn~Y s » {d) " cl ~ {cl ~ cl)u °-° 1 2 ' ( 1 ■ 1 0 ) 

which is a nonlinear least squares optimization problem. 
Specifically, the unknown parameters in equation (1.6) can be 
determined by solving 



AM 

min f = V 

9,c\,c2,a J 

n=0 



y n -^s n {0)-c\-{c2-c\)u n _ a \ 2 . (i.ii) 



To determine the solution of (1.11) , it is useful to 

partition equation (1.11) over a first interval before the 

transition of the square wave and a second interval after the 

transition of the square wave. That is, equation (1.11) 
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becomes 



cr-1 



/-I 



n=0 



N-\ 



+ 



(1.12) 



Equation (1.12) may be minimized over [0 , cl, c2, a) 
jointly using any of the previously mentioned optimization 
methods. However, since cl and c2 are in separate terms of the 
objective function, their estimates also may be solved for 
separately. For example, the estimate for cl may be obtained 
analytically by differentiating the portion of equation (1.12) 
that corresponds to the first interval with respect to cl, 
setting the result equal to zero, and solving for cl . The 
estimate of c2 may be solved by operating upon the portion of 
equation (1.12) that corresponds to the second interval in like 
fashion. 

The estimates of cl and c2 also may be obtained 
qualitatively. For example, the estimate for cl may be 
expressed as a mean of an error between the observation y n and 
the signal prediction s n {6) before the square wave transitions; 
similarly, the estimate for c2 may be expressed as a mean of an 
error between the observation y n and the signal prediction s n {9) 
after the square wave transitions. Hence, the estimates cl of 
cl and c2 of c2 may be expressed as 



i a~\ 



(1.13) 
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and 



c2 



-_L_Y 



(1.14) 



Equations (1.13) and (1.14) then may be substituted back 
into the objective function of equation (1.12), resultingin 
the following expression of the objective function: 



C* 1 

/-I 



»=0 



l ^ 4, l 



+1 



1 



N-a 



W-a ~ 2 



(1.15) 



Now, equation (1.15) is a function of the observation, y n , 
the unknown signal parameters, 0, and the location of the step 
function, a. All of these parameters may be jointly estimated, 
for example, using non- linear optimization techniques as 
described above. 

Nevertheless, it also may be possible to estimate only the 
unknown parameters cl, c2 and a based on expanding and 
rearranging the terms of equation (1.15) to give 



n=0 



-g(a) 



-a 



a n=0 



1 ° 



-(N-a) 



i 



N-a tt 2 



L a m=0 



+ 2aRe 
+ 2(N-a)Re 



N-a 



N-a 



1 W-a i Af-or ^ 



N-a fa 2 
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(1.16) 



where 



a r; 



+ (N-a) 



n=0 



N-a 



(1.17) 



The first term in equation (1.16) is an expression of mean 
square error between the observation y n and the signal 
A 

prediction -^-s tt {0) , while the second term is explicitly provided 

by equation (1.17) . 

All of the other terms of equation (1.16) involve averages 

A 

of the signal prediction -^-s w (<9)and may be approximated as zero 

if the expectation of s n {9) is approximately equal to zero, for 
both before and after the transition of the step function. When 
the expectation of s n {6) may not be approximated as zero, the 
parameters may be estimated using a method that retains these 
terms. For example, the parameters may be estimated by starting 
with a seed value of a that may be used to determine an 
estimate of 0, which, in turn, may be used to produce an 
estimate of a. The method may be iterative and may continue to 
alternate between estimating a and 6 until convergence to a 
desired degree of precision is achieved. 

Nevertheless, for many signals, such as, for example, a GSM 
signal for which the expected value of the underlying binary 
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data stream is zero or approximately zero, it is reasonable to 
assume that the average of the signal prediction s n (6) is equal 
or approximately equal to zero when taken over a sufficiently 
large interval. For example, the signal prediction may be 
expressed as: 



k=0 



(1.18) 



where d k is an original binary data sequence with an expectation 
of zero, and h n is the combined action of the transmit filter, 
the channel filter, and the receive filter. Here the vector of 
unknown parameters, 0, can be taken as the complete data 
sequence d k for all k and the complete filter h n for all n. 
Because the expectation of the binary sequence is zero and the 
binary sequence is independent of the combined filter, then the 
expectation of the signal in equation (1.18) is zero. That is, 

E{d k } = 0 implies = 0. (1-19) 

Hence, approximating as zero the expectation of s n {6) , the 
objective function of (1.16) becomes: 



AM 



n=0 



(1.20) 



and equation (1.20) may be minimized by selecting an a 
maximizes g(a) . That is, 



that 



a « argmaxg(a) = argmaxa 



+ (N-a) 



1 

N-a 



N-a 



(1.21) 
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Fig. 5 illustrates a method 139 for optimizing equation (1.21) 
that may be used to implement the system of Fig. 1. A sum of the 
received data is computed and stored (step 510) , where the sum may 
be expressed as : 



AM 



*>5>. ■ d-22) 

Next, temporary parameters Y p3 (a partial sum of the data) 
and g max are set initially to zero, and temporary parameter a Test 
is set initially equal to one (step 520) . 

Using the parameters of step 520, estimates g m3LXI a, and Y ps 

may be computed iteratively over increasing values of a Tes t while 
a Test is less than or equal to N-l, the number of data samples 
(steps 530) . For example, as shown in Fig. 5, estimating g max/ 
d , and Y may include adding the current data Y„ to the 

F a Tesi-\ 

partial sum of the data Y ps to generate an updated partial sum Y ps 
(step 533) . An updated value for the object function g then may 
be determined as: 

8= ^-\ Y »\ 2+ T^r\ Y *- Y »\ 2 (l - 23) 

u Test iy u Test 

(step 535) . The updated value of g then may be compared to the 
stored value of g max (step 537) , and if updated g is greater than 
9max, then g max may be set equal to updated g as a best current 
guess of the maximum of g, d may be set equal to a Te3tf and Y 

may be set equal to Y ps (step 539) . After updating the values of 
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gmax, d, and Y ps (step 539), a Test may be incremented (step 541) 

and, if a Test is less than or equal to N-l (step 531) , then steps 
53 0 may be repeated. 

The estimation of the parameters accomplished in steps 530 
also may be performed, for example, by decrementing a Test from a 
high value to a low value, or by performing a random selection 
of a Test . Under any of the approaches mentioned, parameters may 
or may not be estimated for each value of a Test . 

Following completion of the iterative process of steps 530, 

the estimated values of g maXf d , and Y ps may be used to correct 

the offset of the data y n (step 550) . For example, using d as 
the estimate of the transition point of the step function, the 
estimate c\ may be expressed using the calculated values as 

£1 = 4^* (1.24) 
a 

while 62 may be expressed as 

£2= JTd^ - f »}- (1 - 25) 

Optionally, where cl and c2 as estimated above are equal or 
approximately equal (indicating that a step function may not be 
present) , then both may be re-estimated as 

cl = c2 = -l(rj. (1.26) 

Thereafter, to correct the offset of the received data y ni the 
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estimated parameters may be used to subtract the step function 
from each data point as follows 

ty.-Sl, 0<n<d 
v = < (1 

y * \y n -c2, «<«<A^. 
Following the correction of the offset, further estimation 
methods may be applied to the residual data {y n minus the step 
function) in order to estimate the remaining unknown signal 
parameters 9. 

Other implementations are within the scope of the 
following claims. 
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